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chapter    1 
INTRODUCTION 

This   report    describes    the    status   of    an    investigation    into 
the   behavior    of    a   composite  structure    in    a    fire    environment. 
The    investigation   arises    from   the   need    to    assess    the    damage 
a    composite   aircraft    structure   undergoes    when   exposed    to    a 
carrier    deck    pool   fire.      Towards    this    end,    the    Naval    Weapons 
Center    (NWC)     in    China    Lake    has   undertaken    a    series   of    exper- 
iments  to    determine   the   effects    of   fires   on    structures     (1). 
As    the  number   of   parameters   eff acting   thermal    behavior    is 
large,    and   as  the   behavior    is    inherently    nonlinear,    an    ex- 
haustive   experimental   study  is    prohibitive    in    cost    and   rime. 
To    assist    their    experimental    program,    NWC    has    funded    the 
NAVAL    POSTGRADUATE    SCHOOL     (NP S)     to   develop    a    mathematical 
model   of    the    problem.      The    NWC   experimental    investigation, 
and    the   NPS    analytical    investigation    are    intended    to    comple- 
ment  one    another. 

The  NPS    effort    seeks    to    quantitatively    predict   the    extent 
of    damage    a   composite    structure    undergoes    when   exposed    to    a 
fire.      This   objective   was    to   be    achieved    through    the    devel- 
opment  of    a   mathematical    model    of   the    problem.      The   computer 
program   resulting    from   the    mathematical    model    would    be    used 
to    predict    the   thermal    behavior    of   composites    for    a   wide 
range   of    system    parameters.      The    parameters    characterizing    a 
system   include    (1)    composite   components    (materials),     (2) 
composite    layup     (lamina    orientation),     (3)     composite   thick- 
ness,    (4)     permeability,     (5)    porosity,     (5)     wind,    and    (7)    fire 
intensity.      The    degree   of    validity   of   the    model    could    be    es- 
tablished   by    comparing    the    results   of    the    model   to   experi- 
mental results.       Moreover,    the    model    would    be    used   to    sug- 
gest   the    most    appropriate    experiments. 

The  overall    problem    of    the    behavior   of    a    composite    in    a 
fire    environment    is   conveniently    partitioned    into    two    inves- 
tigations;   an   investigation   of   thermal    behavior,    and    an 
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investigation  of  how  composing  strength  is  affected  by  it's 
temperature  distribution. 

Mathematical  models  of  both  thermal  behavior  and  strength 
behavior  have  been  developed.   The  computer  codas  resulting 
from  the  mathematical  models  have  been  implemented  on  a  num- 
ber of  cases.   Some  results  and  brief  descriptions  of  the 
thermal  and  strength  analyses  follow. 
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chapter    2 
THERMAL    INVESTIGATION 


The  thermal   investigation   is    concerned    with   the   determi- 
nation of    the   temperature    profile    in    a    coaposite    laminate 
due    to  exposure    to   a    fire    environment. 


2. 1  OBJECTIVES 

The   objectives   of    the   thermal    investigation    are: 

1.  development  of  a  mathematical  model  for  the  determi- 
nation of  behavior  of  a  composite  laminate  in  a  fire 
environment. 

2.  to  shew  how  each  of  the  system  parameters  affect  the 
thermal    behavior. 

3.  to    evaluate   the    extent,    and   rate    of    material    consump- 
tion  due    to   combustion. 

4.  to    extend   the    Semenov   combustion    model    to   composite 
media. 

2.2  DESCRIPTION    OP    THE    MODEL 

The   principal    investigator    for   the    determination    of    ther- 
mal   response    of    a    composite   laminate    in    a    fire    environment, 
until    the    receipt    of    his    Doctoral   Degree    in    June    1932,    has 
been    Lcdr.    C.5.    Vatikiotis    of    the   U.S.    Navy.       Initaily    Lcdr. 
Vatikiotis    (then   a    Lieutenant)    worked    on    the    problem    to    sa- 
tisfy the    research   requirement    towards   an    Sr.aineer 's    Degree 
in    Mechanical   Engineering   from    LJPS.      The    iegree    was    granted 
in    June    1S80.      A   brief    description   of    the    mathematical    model 
for    the   thermal    behavior    follows.      A    more    detailed    descrip- 
tion   is    contained   in    references     (2,3)  . 

We   consider   a   composite    of    graphite    fiDers    imbedded    in    an 
epoxy   matrix.      As    epoxy    combusts    at   a    much    lower    temperature 
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than    graphite,    the   model   assumes    the   epoxy    has    burnt    off   and 
all    that    remains  is  a   porous   graphite   medium.      Figure    1    pre- 
sents   a    diagram    of   the    geometry    of   the    porous   medium.       A 
pressure    differential   across  the    porous   plate   induces    an   air 
flow   through   the   medium.       The    air   flow    provides    two   opposing 
actions;    heat    generation    and    heat   transfer.       The   heat    gener- 
ation  results   from   the   combustion   which    requires   oxygen   to 
be    ever    present.       Heat    is    transferred    thrsugh    the    porous 
medium  by    convective   air    flow.       Behavior    of    the    system    de- 
pends  upon    which   of   these    effects,    heat    generation    or    heat 
transfer,    dominates. 

A    one    dimensional    model    of    a    porous   medium    was    developed 
on    the   basis    of   the    following    fundamental    laws: 

1.  an    energy   balance   on   the    graphite, 

2.  an    energy   balance    on  the    air,    and 

3.  a    mass    balance   on   the    oxygen. 

All  cf  the  heat  transfer  modes,  conduction,  convection,  an  I 
radiation  are  included  in  the  no  del.  In  addition,  the  ener- 
gy balance  equation  for  the  graphite  contains  a  heat  genera- 
tion term  of  Arrhenius  type  arising  from  combustion  of  the 
graphite.  similarly,  the  oxygen  mass  balance  equation  con- 
tains a  term  accounting  for  the  consumptian  of  oxygen  due  to 
combustion . 

In   addition    to   the   energy    and    mass    balance    equations, 
Darcy's    law    for    rhe    governing    equation    for    air    flow   through 
the    porous    graphite    plate    was    combined   wi-h    Bernoulli's 
equation    to   yield   an    equation   for   the    pressure    differential. 
As    very   large  changes    in    temperature    occur    during   the    tran- 
sient  history    of   the   system,    all    thermophysicai    properties 
are    treated   as  temperature    (and    hence   time)     dependent.       The 
combustion    model    was    patterned   after    the    Semenov    model,    and 
the    reaction    was   taken   as    producing   only    oarbon    dioxide. 

The   mathematical    model    of   the    thermal   problem   consists    of 
three   nonlinear,    coupled,    partial   differential    equations    in 
time   and    space.      These   field    equations    were    transformed,    via 
a    GaierJcin    formulation    of   the    finite    element    method,    into    a 
set    cf  3 xn   nonlinear,    coupled   ordinary    differential   equa- 
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Figure  1.  Geometry  of  the  Model 


tions  in  time,  where  n  is  the  number  of  spatial  ncdes  in  the 
finite  element  discretization.   These  equations,  upon  numer- 
ical integration  by  the  Franke  (4)  modified  Gear  method  for 
stiff  systems,  gives  the  graphite  temperature,  the  air  temp- 
erature, and  the  oxygen  concentration  at  each  nodal  point  as 
a  function  of  time.   In  addition  to  these  response  varia- 
bles, a  large  number  of  other  system  parameters  are  calcu- 
lated.  These  include  air  density,  porosity,  permeability, 
Reynold's  number,  pore  velocity,  and  pressure.   In  this  way 
a  full  characterization  of  the  system  is  known  during  the 
transient  history  of  the  problem. 

2.3    GEHEBALIZATIOH  OP  THE  MODEL 

On  the  basis  of  his  NPS  research  activity,  Lcdr.  7ati- 
kiotis  was  accepted  to  the  NPS  Doctoral  program.   He  was 
permitted  to  conduct  the  research  require!  for  the  Doctoral 
Degree  at  his  next  duty  station,  the  Davizl  Taylor  Naval  Ship 
Pesearch  and  Development  Center  (NSRDC) .   During  the  period 
of  June  1980  to  June  1982,  Lcdr.  Vatikiotis  made  significant 
improvements  to  the  thermal  modal,  and  performed  a  number  of 
computer  runs  to  determine  the  affects  of  design  and  envi- 
ronmental parameters  on  system  behavior.   A  brief  descrip- 
tion of  his  work  is  presented  nere.   For  a  acre  detailed  ac- 
count of  the  work,  the  reader  is  referred  to  the  Doctoral 
dissertation  (5)  . 

2-3.1    Modularization  of  the  Computer  Program 
The  computer  program  generated  in  the  Engineer's  Degree 
research  effort  was  put  into  modular  form  in  order  to  more 
easily  accommodate  any  changes  to  the  model.   At  present, 
the  computer  code  consists  of  25  subroutines. 

2.3.2    Transient  Porosity  and  Permeability 

In  the  earlier  model,  porosity  and  permeability  were 
treated  as  uniform  through  the  medium,  and  constant  with  re- 
spect to  time.   The  model  was  updated  to  account  for  an  ini- 
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rial    non    uniform    distribution   of    graphite    filaments.       More- 
over  as   graphite   is   consumed   by    combustion,    both    the    porosi- 
ty   distribution    and    permeability    distribution    vary    with   time 
and    position.      The   model   was    modified    to   account    for   the 
changing    air   density,    and    porosity  by    invoking   the  continui- 
ty   equation.      Combining    the  continuity   equation    with    Darcy's 
law    resulted    in    a    fourth   nonlinear   second    order    field    equa- 
tion.     This   equation    was   solved    numerically    by   the    "shooting 
method". 

2.3.3  Combustion 

In  the  earlier  model  the  product  of  ths  graphite  combus- 
tion was  assumed  to  be  carbon  dioxide  only  (i.e.,  a  simple 
reaction).   The  literature  shows  that  carbon  monoxide  be- 
comes a  byproduct,  along  with  carbon  dioxide,  at  higher  com- 
bustion temperatures.   An  empirical  relation  for  the  propcr- 
tion  of  CO  and  C02  as  a  function  of  temperature  was 
incorporated  into  the  model. 

2.3.4  Arrhenius    Expression 

As   Prank-Kamenetskii     (6)    suggested    that    carbon    coabustion 
is    net  a    first   order    reaction,    the   model    was    generalized    to 
accommodate    fractional    order   reactions.      This    modification 
was    accomplished   by   raising  the    oxygen    concentration    terni    in 
the    reaction   rate    expression    to    a    fractional    power.       A    ser- 
ies   cf  analyses    were    performed  with    fractional    exponents    1/2, 
2/3,    and    3/u.      The   analytical    results   changed    markedly    with 
exponent.       Hence,    the   accuracy   of   any    analytical    model    de- 
pends  on    the    accuracy   of   the    experimentally    derived   exponent. 


2. 3.5        Transient    Plate   Thickness 

Accounting    for    the   change   in    porosity    iue    to    graphite 
consumption    also   permitted   the   determination    of    the   laminate 
thickness    with   time. 
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2-3-6    Porous  Model 

Originally,  the  porous  medium  was  modeled  as  cylindrical 
fibers.   The  model  has  been  modified  so  that  it  can  treat 
either  spherical  particles  or  cylindrical  fibers.   This  op- 
tion permits  the  user  to  consider  the  problem  of  combustion 
of  graphite  particles. 

2-3.7    Heat  Transfer  Coefficient 

Alternate  models  for  internal  heat  transfer  coefficient, 
surface  area  per  unit  volume,  and  average  pore  diameter  are 
available  for  selection  cf  the  investigator.   in  addition, 
the  pore  diameter  is  based  on  the  hydraulic  diameter. 

2-3-8    Surface  Recession  Problem 

In  the  case  of  combustion,  the  reaction  invariably   moves 
to  the  source  of  the  oxygen,  namely  the  air  inlet  surface. 
This  results  from  the  fact  that  in  a  relatively  short  time, 
all  cxygen  in  the  interior  is  consumed.   Dnce  all  of  the  in- 
terior oxygen  has  been  consumed,  all  subsequent  combustion 
takes  place  at  the  air  inlet  surface  where  there  is  always  a 
continuous  supply  of  cxygen.   At  this  stage  in  the  transient 
history,  the  problem  degenerates  intc  a  "surface  recession" 
or cblem.   Since  there  is  no  oxygen  except  at  the  inlet  sur- 
face, the  oxygen  balance  equations  must  be  removed  from  the 
analysis,  otherwise  severe  numerical  difficulties  arise. 
The  removal  cf  these  equations  from  the  analysis  is  accom- 
plished automatically  in  the  computer  program,  so  that  the 
user  need  not  be  concerned  with  the  transformation  to  a  sur- 
face recession  problem. 


2.3.9    Initial  Conditions 

The  arbitrariness  of  selecting  initial  conditions  was 
overcome  by  incorporating  into  the  analysis,  a  subroutine 
which  generates  a  transient  thermal  environment.   If  the 
user  selects  this  mode  for  the  initial  conditions,  a  heat 


flux  and  it's  duration  tint 
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specified  by  the  use; 
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2.3.10  Boundary   Conditions 

The  model    was   given   more  flexibility   for    application   to   a 
wider   class   of    problems    by    allowing   a    choice    of    boundary 
conditions.      The   user   may    select    any    of   the    following    three 
boundary    conditions: 

1.  insulated   surfaces    and   DankwertJs    B.C.'s    for   the   air 
temperature   and   oxygen   concentration. 

2.  radiation  at   surfaces   and    Dankwert's    B.C.'s    for    the 
air  temperature   and    oxygen  concentration. 

3.  convective   and   radiative    heat    transfer    at    both    sur- 
faces,   air    temperature    at    inlet    surface   equal    to    am- 
bient  temperature,    oxygen    concentration    at   the    inlet 
surface    equal    to    ambient    concentration,    and   Dank- 
wert's  B.C.'s    at    the   exit    surfaces. 

2.3.11  Numerical   Instabilities 

During  certain  times  in  the  transient  history,  numerical 
instabilities  were  encountered.   This  is  particularly  true 
when  the  oxygen  concentration  is  going  to  zero  in  the  inte- 
rior of  the  medium.   The  numerical  instability  appears  in 
the  results  as  an  oxygen  concentration  oscillating  about 
zero.   A  subroutine  was  added  to  the  program  which  dampens 
the  oscillation  significantly. 

2.3.12  Improved  Accur acy 

In  the  original  formulation  of  the  problem,  the  nonlinear 
heat  generation  term  vas  treated  as  a  "forcing"  term.   A 
significant  reduction  in  computational  effort  resulted  when 
this  nonlinear  term  was  treated  as  a  part  of  the  field  oper- 
ator.  This  improvement  in  computational  effort  is  attribut- 
ed to  the  fact  that  as  an  operator,  the  generation  term  en- 
ters the  analysis  as  part  of  the  Jacobian  matrix  of  the 
system.   As  the  Jacobian  matrix  plays  a  prominent  role  in 
the  numerical  integration,  it  is  no  wonder  that  a  signifi- 
cant improvement  occurs. 
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2.3.13  Updating  Theriop.hysical  Properties 

In  the  original  formula tion,  all  temperature  dependent 
properties  were  updated  outside  of  the  time  integration 
loop,  at  the  arbitrary  discretion  of  the  iser.   The  formula- 
tion was  modified  so  that  all  temperature  dependent  proper- 
ties are  updated  continuously  and  automatically  within  the 
integration  loop.   The  result  is  improved  accuracy  ar.d  a  re- 
duction of  computational  effort. 

2.3.14  Thermal   Conductivity 

The  expression    for   thermal    conductivity    of   the    porous 
medium  was    changed   to   one    that    axtended   the    range   of    porosi- 
ties   which    could   be   treated.      The   magnitules   of    the   internal 
parameters,    such   as    pore    velocity,    etc.,    are    checked   to   ver- 
ify   that    the   empirical    expressions   used    to    calculate    these 
parameters    remain    valid. 

All   of    the   modifications   briefly   described    above   had   the 
effect  of     (1)    extending    the   raoiel   to    a   more    general   state, 
(2)     provide   the    analysis    with    more   accuracy,    and    (3)    reduce 
the    computational    effort    required    in    the    solution    of    a    prob- 
lem.     In    the   next   section    some   of   the    results    obtained    by 
Lcdr.      Vatikiotis    are   briefly    described. 

2.4    DISSERTATION  RESULTS 

A  number  of  computer  analyses  were  performed  in  order  to  ob- 
tain results  which  would  provide  some  understanding  of  the 
effects  of  each  of  the  system  parameters  on  system  behavior. 
In  order  tc  establish  the  validity  of  the  mathematical  mod- 
el, a  comparison  of  results  from  computer  analyses  with  ex- 
perimental results  was  made  whenever  possible.  In  such 
cases  there  was  fair  agreement  between  analytical  and 
experimental  results.  There  are  three  reasons  why  excel- 
lent agreement  between  experimental  and  analytical  results 
should  not  be  expected.  First,  combustion  experiments  by 
their  very  nature  are  difficult  to  perform  and  control. 
Thus  there  is  a  large  uncertainty  in  the  reported  results. 
The  relatively   small  number   of  combustion   experiments  at- 
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tests  to  their  difficulty.  Secondly,  there  is  small  likeli- 
hood that  the  values  of  many  system  parameters  is  known  to 
any  high  degree  of  certainty.  Thus,  for  example,  the  expo- 
nent  en  oxygen  concentration  ia  the  Arrhenius  expression  has 
not  been  accurately  determined;  and  yet  a  small  change  in 
the  exponent  yields  significant  changes  in  the  quantitative 
results.  Thirdly,  there  is  the  inherent  uncertainty  associ- 
ated with  every  complex  mathematical  model.  Any  number  of 
simplifying  assumptions  are  made;  and  of  course  there  are 
choices  in  the  selection  of  competing  relations  for  calcu- 
lating properties.  Ail  of  these  considerations  suggest  that 
the  results  be  regarded  as  descriptive  rather  than  precise. 
As  ncne  of  the  analysis  results  contradicted  experimental 
results,  there  is  confidence  that  the  mathematical  model  is 
a  good  one.  Figure  2  presents  a  comparison  of  the  mathemat- 
ical model  and  the  experimental  results  of  Kolodstev  (7) . 
The    data    for    the   analysis    is   given   in    Table    1. 

A    series   of   computer    runs    were    performed    to    isolate    the 
effects    of    a    single    system    parameter    on    system    behavior.       A 
brief   description    of    each    of   the    studies    follows.      In    each 
case,    the    boundary    conditions    with   insulated    surfaces    was 
used.      7'ae    numerical    results   given    in    this    report    should    be 
considered   approximate    rather    than    precise. 

2.4.1        Effects   of   Fire   Fnyironm  ent 

In   this    study,    the   effect    of    fire    intensity   on    a    compos- 
ite   plate      was   investigated   by   subjecting    a    plate    to    two 
heat   flux   conditions.      The    first    heat    flux    condition    results 
in    combustion,    and   the    second    heat   flux    condition   results    in 
extinguishment.      Table    2    presents   the   data    describing    the 
composite   of   this    study. 

When    a    heat    flux   of    30,000    Stu/ft-hr    is    applied   to    the 
inlet   surface   for    15    seconds,    and    is    then    removed,    the    plate 
undergoes    combustion.       Figure    3    shows    the    graphite    and    air 
temperatures    as    functions    of    space   at    selected   times.       Fig- 
ure   4    plots   oxygen    concentration      versus    space    at    selected 
times,    and   Figure    5    shows   the    thickness    of    the    plate    as    a 
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TABLE  1   Geometry  of  porous  medium  and  ambient  conditions 

■for  simulating  Kolodstev's  C  7  1    experiment. 


•  Particle  Shape:  spherical 

•  Particle  Diameter  (  d  )  :  .126  in 
Unit  Cell  Thickness  <D>:  .126  in 

•  Spatial  Thickness  o-f  Porous  Medium  (L  ):  7.5  in 

Porosity  ( p  )  . 476 

-6     2 

Permeability  (  m )  .103x10    -ft 

Bulk  Thermal  Conductivity  o-f  Carbon  Ckc):  36.0  Btu/ft-hr — F 

Bulk  Specific  Heat  o-f  Carbon  (  Cc )  :  .231  Btu/lbm-F 

Bulk  Density  of  Carbon  (p  ):  70.3  lbm/'ft 

Thermal  Emissivity  of  Particles  (€):  .9 

Ambient  Temperature  (  7W)  :  SO.O  deg-F 

•  Ambient  Pressure  (  F^>  :  14.7  psi 
•Pressure  Differential  (aP) :  -.29  psi 

•  Ambient  Oxygen  Concentration  (CD.):  .0172  lbm/ft 


•Known  parameter  for  Kolodstev's  experiment. 
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function    of  time.      Figure    4   shows    that   the    graphite   and   air 
tem  peratures    differ   only   at   the    air   inlet    surface    (X/L=0.) 
until   about    t=360    sec.      Thereafter   the   air    temperature    is 
sign  if icantly   less   than    the  graphite    temperature    for    a    con- 
siderable   range    of   X/L.      This    can   be    misleading    because   the 
thickness    L   has    changed    significantly   during    the    evolution. 
Figure  5    shews   that   the    plate   thickness    at    5  17    seconds    is 
only    10    percent    of    it's    original    thickness.       Thus   the    pene- 
tration   of    "cooler"    air    is    .0025    inches.       At    20    seconds,    the 
plate   thickness    is    .25    inches    and  the   cooler    air    penetrates 
•"The    plate    to    X/L=0.1,    or    .025    inches.      Thus    when   looking    at 
Figure   3,    one    must   keep    in    mind    that    L,    the    thickness    of   the 
plate,    is    changing    due    to   combustion. 

An   important    observation   which   can   be    oiade    for   this    com- 
bustion  case    comes   from   Figure   '4  .      Figure    4    shows   that    the 
oxygen   in    the   interior   of    the    plate   is   very    quickly   con- 
sumed.     All   cases    which    have   gone    to    combustion    have    shown 
this   trait,   the    interior   oxygen    is   quickly    depleted-       At 
time    equal   to   20    seconds,    the    problem   is   about    to    become    a 
"surface    recession"    problem,    in    which    all    reaction    takes 
place   at    the    air    inlet    surface.       Figure    5    shows    that    the 
rate   of   material   consumption   becomes    ever    increasingly 
greater    with   time    after    surface    recession    begins.      For    this 
particular    case   we   note    that   the    temperature    increased    from 
a    maximum    of    900    degrees   Fahrenheit    when    the    heat    flux    was 
removed   to    a   temperature    slightly   less   than    3000    degrees 
Fahrenheit    at   5  1"?    seconds. 

When    a    heat    flux   of    30,000    3tu/ft-hr    is    applied    to    the 
inlet    surface   for    14    seconds   and    then    removed,    the    plate 
goes    to    extinction    in    less    than    a    minute.       Figure    6    shows 
the    manner    in   which  the    temperature   distribution    goes    to    the 
uniform   ambient    temperature.      We    note   that    the    peak   tempera- 
ture   while    continuously    decreasing   in    magnitude    with    time, 
moves    from    the    inlet    surface   to    the   exit    surface.      Figure    7 
shows   that    as   the    plate    cools,    the   oxygen    concentration    in- 
creases   until   it   reaches    the    uniform    ambient    condition. 
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TABLE  2 Geometry  o-f  porous  medium  and  ambient  conditions 

■for  the  example  problems  in  Section  2.4.1 


Particle  Shape:  spherical 

Particle  Diameter  (d):  .005  in 

Unit  Cell  Thickness  (D>:  .005  in 

Spatial  Thickness  o-f  Porous  Medium  (  L>:  .25  in 

Porosity  (  p  )  . 476 

-9    2 

Permeability  <  m  )  .162x10^  -ftc 

Bulk  Thermal  Conductivity  o-f  Carbon  (  kc)  :  86.0  3tu/-ft-hr — F 

Bulk  Speci-fic  Heat  o-f  Carbon  (  Cc  )  :  .231  Btu/lbm-F 

Bulk  Density  o-f  Carbon  C& )  :  70.3  lbm/-ft 

Thermal  Emissivity  o-f  Particles  (€  ):  .9 

Ambient  Temperature  (T  >:  30. 0  deq-F 

CO 

Ambient  Pressure  (  P^)  :  14.7  psi 

Pressure  Di  -f  -f  erenti  al  (Ap)  :  -.35  psi 

3 

Ambient  Oxygen  Concentration  ( <D)  :  .0172  lbm/ft 
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The  point  to  be  made  here  is  that  a  slight  change  in  the 
duration  of  heat  flux  can  determine  whether  damaging  combus- 
tion or  safe  extinction  occurs.   Here  we  have  seen  that 
30,000  3tu/ft-hr  for  15  seconds  results  in  a  combustion 
which  completely  consumed  the  plate.   The  same  plate  goes  to 
extinction  if  the  heat  flux  is  terminated  after  14  seconds. 

2. 4.2   Effects  of  Permeability 

Permeability,  a  measure  of  flew  conductance,  is  one  of 
the  composite  design  parameters  which  most  affects  thermal 
behavior.   In  order  to  isolate  the  effect  of  permeability  on 
behavior,  the  geometry  of  the  filament  diameter  and  the 
width  of  the  unit  cell  were  varied  in  such  a  way  as  to  main- 
tain a  constant  porosity  for  aii  analyses  in  this  study.   In 
this  way,  the  quantity  of  air  around  a  filament  at  any  in- 
stant is  fixed,  but  the  rate  of  air  passing  a  filament  in- 
creases with  increasing  permeability.   Table  3  lists  tne 
data  for  a  number  of  analyses  which  were  used  in  the  perme- 
ability study.   Figure  3  presents  temperature  versus  perme- 
ability for  three  values  of  pressure  differential.   The 
temperature  here  is  the  temperature  above  which  the  plate 
will  experience  combustion,  and  below  which  the  plate  will 
extinguish.   The  "ignition"  temperature  for  each  case  is  ob- 
tained by  trial  and  error  iteration  on  a  uniform  initial 
temperature.   After  a  number  of  runs,  we  obtain  two  nearly 
equal  uniform  initial  temperatures,  one  leading  to  combus- 
tion and  the  other  extinguishment.   Their  average  is  taken 
as  the  "ignition"  temperature. 

Figure  8  shows  that  for  a  given  permeability,  the  temper- 
ature at  which  combustion  occurs  increases  with  increasing 
pressure  differential.   An  increase  in  the  pressure  differ- 
ential causes  an  increase  in  the  air  flow  through  the  porous 
medium.   Here  an  increase  in  air  flow  results  in  "blowout" 
extinguishment.   That  is,  an  increase  in  air  flow  results  in 
an  increase  in  the  ignition  temperature.   We  also  note  that 
the  rate  of  change  of  ignition  temperature  decreases  with 
increasing  permeability. 
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TABLE  3   Geometry  o-f  porous  medium  and  ambient  conditions 

■for  the  permeability  analysis  in  Section  2.4.2 


Particle  Shape:  spherical 

Particle  Diameter  (d  >:  (various) 

Unit  Cell  Thickness  (D>!  (various) 

Spatial  Thickness  o-f  Porous  Medium  (  L  >  *       1.0  in 

Porosity  ( p )  .476 

Permeability  ' m  )  (various) 
Bulk  Thermal  Conductivity  o-f  Carbon  (  kc)  :     36.0  Btu/ft-hr-F 

Bulk  Speci-fic  Heat  o-f  Carbon  (Q):  .231  Btu/lbm-F 

Bulk  Density  o-f  Carbon  (pe  )  :  70.3  lbm/-ft^ 

Thermal  Emissivity  o-f  Particles  (€  >:            .9 

Ambient  Temperature  (  T^) :  SO.O  deg-F 

Ambient  Pressure  (p):  14.7  psi 

CD  ' 

Pressure    Di-f-fersnti  al     (&P)  :  (various) 

Ambient    Oxygen    Concentration     ((t^):  .0172    lbm/ft 
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2.4.3  Effects  of  Porosity. 

The  affect  of  porosity,  the  measure  of  void  volume  per 
unit  volume  of  space,  on  thermal  behavior  is  shown  in  Figure 
9.   The  data  for  the  porosity  study  is  given  in  Table  4. 
Figure  9  shows  that  for  a  fixed  porosity,  the  ignition  tem- 
perature increases  with  increasing  permeability.   Again  this 
shows  -chat  increasing  the  air  flow  will  cause  blowout  at  a 
fixed  temperature,  or  stating  it  a  little  differently,  in- 
creasing the  air  flow  through  a  plate  makes  it  less  likely 
to  combust. 

However,  qiven  a  fixed  permeability,  the  ignition  temper- 
ature decreases  with  increasing  porosity.   The  reason  for 
this  is  that  increasing  the  porosity  while  fixing  all  other 
parameters  results  in  decreasing  air  flow  through  the  medi- 
um.  We  note  that  the  rate  of  change  of  ignition  temperature 
diminishes  with  porosity.   That  is,  the  decrease  of  ignition 
temperature  with  porosity  becomes  less  with  increasing  poro- 
sity. 

2.4.4  Effects  of  Thickness 

Figure  10  shows  the  effect  of  plate  thickness  on  thermal 
behavior.   Table  5  is  the  data  for  the  thickness  study.   The 
most  important  feature  of  Figure  10  is  the  very  large  de- 
crease in  ignition  temperature  with  thickness,  especially 
for  the  smaller  thicknesses  fr;m  0.25  inches  to  2.0  inches. 
An  expiaination  for  this  oehavior  is  the  following.   If  ail 
parameters  except  thickness  are  fixed,  then  increasing  the 
thickness  while  keeping  the  pressure  differential  constant 
results  in  decreasing  the  pressure  gradient.   The  decrease 
in  pressure  gradient  causes  a  ie crease  in  the  pcre  velocity 
of  the  air  through  the  plate.   Thus,  increasing  thickness 
results  in  a  decrease  of  air  flow.   With  less  air  flow,  con- 
vective  heat  transfer  diminishes  and  there  is  i<=ss  opportu- 
nity for  blowout.   Hence  an  increase  in  thickness  should  re- 
sult in  a  lower  ignition  temperature. 

Another  expiaination  for  iecreasing  ignition  temperature 
with  increasing  thickness  is  ae  follows.   The  above  sentence 
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TABLE   4   Geometry  o-f  porous  medium  and  ambient  conditions 

■for  the  porosity  analysis  in  Section  2.4.3 


Particle  Shape:  spherical 

Particle  Diameter  (d  ):  (various) 

Unit  Cell  Thickness  (D  >s  (various) 

Spatial  Thickness  o-f  Porous  Medium  (L  ):       1.0  in 

Porosity  ( p )  (various) 

Permeability  (  nn  )  (various) 
Bulk  Thermal  Conductivity  o-f  Carbon  (  kc)  :     36.0  Btu/-ft-hr — F 

Bulk  Speci-fic  Heat  o-f  Carbon  (Cc):  .231  Btu/lbm-F 

Bulk  Density  o-f  Carbon  <pc>:  70.3  lbm/-ft 

Thermal  Emissivity  o-f  Particles  (  6  )  :           .9 

Ambient  Temperature  ( T— )  :  SO.O  deg-F 

Ambient  Pressure  (  P)  :  14.7  psi 

CD  * 

Pressure  Di  -f-f  erential  (AP)  :  -.5  psi 

Ambient  Oxygen  Concentration  (d)  > s  .0172  lbm/ft^ 
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TABLE  5   Geometry  of  porous  medium  and  ambient  condition* 

for  the  thickness  analysis  in  Section  2.4.4 


Particle  Shape:  spherical 

Particle  Diameter  (d  ):  (various) 

Unit  Cell  Thickness  (Q):  (various) 

Spatial  Thickness  of  Porous  Medium  (  L  ):  (various) 

Porosity  (  p  )  .476 

Permeability  (m)  (various) 

Bulk  Thermal  Conductivity  o-f  Carbon  (  kc)  :     86.0  Btu/ft-hr — F 

Bulk  Specific  Heat  o-f  Carbon  (  Ce >  :               .231  Btu/lbm-F 

Bulk  Density  o-f  Carbon  CD):  70.3  lbm/ft^ 

Thermal  Emissivity  o-f  Particles  (€  ):  .9 

Ambient  Temperature  (  T^)  :  SO.O  deg-F 

Ambient  Pressure  (  f^) :  14.7  psi 

Pressure  Differential  (Ap) :  -.5  psi 

Ambient  Oxygen  Concentration  (<£>J>:                                    .0172  lbm/ft 
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can  be  turned  around  to  "decreasing  the  thickness  will  re- 
sult in  increasing  tha  ignition  tsmpsraturs" .   Now  the  thin- 
ner the  plate,  the  less  distance  the  air  traverses  in  pass- 
ing through  the  plate.   Since  the  air  temperature  increases 
as  it  passes  through  the  plate,  thin  places  have  a  larger 
fraction  of  their  region  with  cool  air  than  do  thick  plates. 
As  thin  plates  are  more  effectively  cooled  by  the  air  flow, 
higher  ignition  temperatures  are  needed  for  combustion. 

Tha  three  curves  again  show  t  hat.  increasing  the  pore  ve- 
locity by  increasing  the  permeability,  leads  to  blowout,  and 
hence  greater  air  flow  (or,  egui vaiently ,  permeability) 
leads  to  higher  ignition  temperatures. 

Another  striking  feature  of  Figure  10  is  the  rapid  de- 
crease in  the  rate  of  change  of  ignition  temperature  with 
thickness.   A  large  change  in  the  ignition  -emperature  can 
be  effected  by  a  small  change  in  the  thickness  when  the 
plate  is  thin.   On  the  other  hand,  very  little  change  in  the 
ignition  temperature  can  be  effected  if  the  plate  is  very 
thick. 

2. 4 . 5   Effects  of  Pressure  Differential 

Each  of  the  previous  studies  sought  to  ieterrair.e  the  ef- 
fect of  a  design  parameter  on  system  behavior.   These  param- 
eters are,  generally,  under  the  control  of  the  structural 
engineer,  and  thus  can  be  selected  to  minimize  combustion. 
In  contrast,  there  are  external  parameters  such  as  pressure 
differential  (possibly  arising  from  wind  tver  one  surface) 
which  are  difficult,  if  at  all  possible,  to  control. 

Table  6  presents  the  data  for  the  pressure  differential 
study,  and  Figure  11  presents  the  results  of  the  study.   The 
results  show  that  the  ignition  tenperatur?  increases  with 
increasing  pressure  differential,  and  that  the  rate  of  in- 
crease decreases  with  pressure  differential.   As  in  the  pre- 
vious studies,  it  is  again  eviient  that  an  increase  in 
permeability  results  in  higher  ignition  temperatures.   The 
-hree  curves  in  Figure  11  show  that  the  effect  of  pressure 
differential  is  much  greater  for  laraer  values  of  permeabil- 
ity. 
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2.4.6        Effects    of    Boundary.  Conditions 

All  of   the  the   previous   parameter   studies -were   performed 
using   insulated    boundary   conditions   at   the    inlet    and    exit 
surfaces.      In  this    study,    the    effect    of   heat    transfer    at    the 
boundaries    was   investigated.      Figure    12    compares   the   behav- 
ior   of  plates   with   insulated   boundaries    (curve    3)    to    plates 
with    heat    transfer   at   the    boundaries     (curve    A).       Although 
the    ignition   temperatures    for    the   plates    with    heat    transfer 
are    always    greater    than    the  ignition    temperatures    for    the 
insulated    plates,    two    regions    are   distinguishable.      In    the 
regicn   with   permeability   greater    than    0.2x10-3    in.    ,    the    ig- 
nition temperatures    for    insulated    plates    are    but    slightly 
less   than    for   "he    heat    transfer    plates.      In    the    region    with 
permeability    less   than    0.2x10-8    in.    ,    there    is    a    far    more 
pronounced    difference   between    irisulated   plates   and   plates 
with    heat    transfer.      The    results    show   that    as    permeability 
decreases,    the   ignition    temperature    decreases    in    both    cases, 
but    the   rate   cf    decrease   is   larger   for    the    insulated   case, 
and    thus    the   difference    between    ignition    temperatures    be- 
comes  increasingly   greater.      Finally,    at    some   small   value    of 
permeability,    the    ignition    temperature   for    plates    with    sur- 
face   heat    transfer    reverses   and    begins    increasing   rapidly 
with    decreasing    permeability. 

An   increase    in    the    ignition   temperature    with    decreasing 
air    flow     (permeability)     for   plates   with   heat    transfer    bound- 
ary   conditions   could    be    explained   if    the    heat    generation    by 
the    exothermal   reaction    ,    and    not   the   heat    transfer    mechan- 
ism,   was    the    dominating    mechanism.      In   this    case    the    reduc- 
tion   of   oxygen    resulting   from    5.    decrease    in    permeability 
would   extinguish   the    combustion    by   "choking".      In   the    chok- 
ing   region,    ignition    temperatures   should    increase   with    de- 
creasing   permeability. 

Apparently,    porous    plates    with   heat   transfer   boundary 
conditions   have    a    minimum    ignition   temperature   at    a    particu- 
lar   value    of   air    flow,    say    Q.      Increasing    the    air    flow    above 
U    results    in   combustion    extinguishment   by    blowout   and    hence 
greater    ignition   temperatures.       This    is   the    case    when    heat 
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TABLE  6    Geometry  of  porous  medium  and  ambient  conditions 

■for  the  pressure  di  f  -f  erenti  al  analysis  in 
Section  2.4,5 


Particle  Shape:  spherical 

Particle  Diameter  (d  ):  (various) 

Unit  Cell  Thickness  (D**  (various) 

Spatial  Thickness  of  Porous  Medium  (  L ) :  1-0  in 

Porosity  (p)  .476 

Permeability  (  rn  )  (various) 

Bulk  Thermal  Conductivity  o-f  Carbon  (  kc)  :     36.0  Btu/-ft-hr-F 

Bulk  Specific  Heat  o-f  Carbon  (  Cc  >  :              .231  Btu/lbm-F 

Bulk  Density  o-f  Carbon  (p<.):  70.3  lbm/ft^ 

Thermal  Emissivity  o-f  Particles  (€  ):  .9 

Ambient  Temoerature  (T  ):  80.0  deg-F 

CO 

Ambient  Pressure  (  P  )  :  14.7  psi 

CD  ~ 

Pressure  Differential  (AP)  :  (various) 

Ambient  Oxygen  Concentration  ((±>  )  :               .0172  lbm/ft^ 
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transfer  is  the  dominant  mechanism.   On  the  other  hand,  de- 
creasing the  air  flow  below  J  results  in  combustion  extin- 
guishment by  choking  and  hence  greater  ignition  tempera- 
tures.  This  is  the  case  when  the  chemical  reaction  is  the 
dominant  mechanism. 

2.4.7   Effects  of  Pore  Velocity. 

Consideration  of  the  results  cf  the  previous  studies  on 
the  effects  cf  individual  parameters  on  system  behavior  sug- 
gests that  the  most  important  parameter  affecting  behavior 
is  the  pore  velocity.   With  this  in  mind,  the  results  of  all 
computer  analyses  (47  problems)  are  shown  in  Figure  13.   The 
curve  through  these  data  points  is  a  fair  fit,  and  thus 
there  is  reason  to  believe  that  the  ignition  temperature  may 
be  obtained,  with  fair  accuracy,  by  the  equation  of  the 
curve. 

This  concludes  the  summary  of  the  work  done  by  Lcdr.  7a- 
tikictis  at  NSRDC  for  his  Doctor  of  Engineering  Degree.   The 
dissertation,  a  more  general  treatment  of  the  combustion 
problem,  presents  a  more  detailed  discussion  of  the  model, 
numerical  aspects  of  the  formulation,  the  results,  and  con- 
clusions.  In  addition,  there  is  a  discussion  of  a  general- 
ised Semen ov  model  for  combustion  of  porous  media,  as  well 
as  an  extensive  list  of  references. 

2.5    NPS  ACTIVITY 

During  the  past  year,  about  half  of  the  funded  research 
activity  at  N?S  was  devoted  to  the  thermal  investigation, 
and  the  other  half  was  given  to  the  strength  investigation. 
In  this  section,  a  summary  of  the  effort  of  the  thermal  in- 
vestigation at  N?s  will  be  described. 

2. 5.  1    Restart  Cagability 

One  of  the  input  parameters  in  the  thermal  analysis  pro- 
gram (CARBGT)  is  the  time  that  the  user  would  like  the  anal- 
ysis to  stop.   It  is  difficult  to  select  a  reasonable  value 
for  this  parameter.   Selecting  too  large  a  time  is  wasteful 
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of  computer  effort,  while  too  small  a  tins  could  require  re- 
doing the  problem.   To  remedy  this  situation,  a  subroutine 
was  added  to  the  program  to  permit  restart.   The  restart 
feature  permits  the  program  to  start  from  the  results  of  a 
previous  analysis.   All  of  the  input  data  necessary  for  the 
restart  analysis  is  obtained  from  the  final  output  of  the 
previous  analysis.   Restart  can  be  invoked  in  the  analysis 
of  a  problem  as  many  times  as  iesirable.   In  this  way,  a 
complete  analysis  of  a  problem  may  consist  of  several  time 
contiguous  analyses. 

2.5.2    Effects  of  Initial  Conditions 

There  are  two  ways  to  input  initial  conditions.   A  way  to 
obtain  reasonable  initial  conditions  on  temperature  is  to 
specify  a  heat  flux  and  a  duration  time,  as  describe!  in 
section  2.3.9.   On  the  other  hand,  if  the  user  is  interested 
in  the  effects  of  the  initial  temperature  on  thermal  behav- 
ior, he  may  input  any  initial  temperature  distribution. 
This  was  the  method  used  in  this  investigation.   Here  we  ob- 
tained the  effects  of  linear  initial  temperature  distribu- 
tions on  ignition  temperature  in  the  following  way.   The 
slope  cf  the  linear  temperature   and  the  average  temperature 
was  specified.   The  same  problem  was  run  as  many  times  as 
needed  tc  obtain  the  average  temperature  below  which  extinc- 
tion occurs  and  above  which  combustion  occurs.   Figure  14 
presents  the  results  of  this  study.   Two  features  of  this 
fiaure  are  (1)  the  maximum  ignition  temperature  occurs  fcr  a 
negative  temperature  slope,  and  (2)  the  ignition  tempera- 
tures are  greater  for  the  negative  temperature  slopes.   Fig- 
ure 14  shows  that  combustion  is  less  likely  to  occur  for 
plates  with  the  higher  temperature  on  the  X/L=0.  surface. 
This  is  a  reasonable  expectation  since  the  cooling  air  is 
most  effective  at  it's  entry  into  the  medium. 
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2.5.3   Graphics  Capability 

Each  computer  analysis  produces  a  vast  quantity  of  numer- 
ical output.   Thus,  at  preselected  instances  of  time  and  at 
each  spatial  nodal  point,  the  program  reports:  (a)  graphite 
temperature,  (b)  air  temperature,  (c)  oxygen  concentration, 
(d)  reaction  rate,  (e)  filament  diameter  (f)  porosity,  (g) 
permeability,  (h)  pressure,  (i)  pressure  gradient  (j)  pore 
velocity  (k)  Reynolds  Number,  (1)  convection  coefficient, 
and  (m)  area  to  volume  ratio.   In  addition,  plate  thickness 
is  reported  at  each  time.   This  data  must  be  assimilated  by 
the  program  user  if  full  benefit  of  the  analysis  is  to  be 
obtained.   As  a  typical  analysis  consists  of  30  to  60  time 
steps  and  25  to  40  spatial  nodal  points,  a  typical  output 
contains  anywhere  from  10,000  tO  30,000  numbers.   It  is  dif- 
ficult, if  at  all  possible,  to  assimilate  this  quantity  of 
numerical  data,  and  therefore  graphical  output  is  more  than 
just  desirable,  it  is  necessary.   To  this  end,  graphics  ca- 
pability was  incorporated  into  the  program.   A  variety  cf 
graphics  options  are  now  available.   The  user  may  select  any 
combination  cf  the  following  output  modes: 

1.  numerical  results  cf  ail  system  parameters, 

2.  curves  of  graphite  temperature  and  oxygen  concentra- 
tion as  functions  of  space  at  preselected  time  in- 
stances, 

3.  three  dimensional  surfaces  of  graphite  temperature, 
oxygen  concentration,  and  reaction  rate  as  functions 
of  time  and  space, 

4.  contour  graphs  of  graphite  temperature,  oxygen  con- 
centration, and  reaction  rate. 

Any  combination  cf  these  options  is   available  at  a  CRT  ter- 
minal or  as  hard  copy  output. 

Two  programs  exist  for  the  surface  and  contour  graphics. 
GRA?3D  produces  surfaces  and/cr  contours  from  the  actual 
numerical  output  of  the  analysis  program,  while  GRAF33 
achieves  smoother  graphics  by  passing  the  output  from  CAR3GT 
through  a  cubic  spline  interpolation  subroutine  (IMSL  sub- 
routine IBCIEU)  . 
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After    inspecting   the    graphical    output    from    GRA73D    or 
GRAF3E,    it    may    be   desirable   to   have   a    mors    detailed   surface 
or    ccntour   of   a    local   subregion    of   the   time-space   domain. 
This    may    be   done    by   recalling    either    graphics   program,    and 
specifying   the    upper    and    lower    limits    of    ths    time   and    space 
subregion    of   interest.       Figures    15   through    20    a:c    sample 
output  of    surfaces   and   contours    from    GRAF3D.       Figures    21 
through    26    ,    are    samples   of  subregion    graphics   of    the    same 
analysis    as   Figures    15    through   20. 

This   concludes   the   discussion    of    the   research   on   the 
thermal    investigation.       A    description    of    the    strength    inves- 
tigation   follows. 
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Figure    15.    Temperature   Surface    from  GRAF3B 
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Figure  16.  Concentration  Surface  from  GRAF3D 
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Figure  17.  Reaction  Rate  Surface  from  GRAF3D 
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Figure    18 .    Temperature   Contours    from   GRAF 3D 
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Figure    19.    Concentration   Contours    from   GRAF3D 


44 


1602 

TIME   (SEC) 


2670 


Figure  20.  Reaction  Rate  Contours  from  GRAF3D 
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Figure  21.  Subregion  Surface  of  Figure  15 


46 


Figure  22.  Subregion  surface  of  Figure  16 
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Figure  23.  Subregion  Surface  of  Figure  17 
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Figure    24.    Subregion   Contour   of    Figure    18 


49 


638  859 

TIME   (SEC) 


1301 


Figure  25.  Subregion  Contour  of  Figure  19 
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Figure    26.    Subregion   Contour   of    Figure    20 
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chapter    3 
STRENGTH   INVESTIGATION 

When    a    composite    laminate    is    subjected    to    a   thermal    dis- 
tribution,   thermal   stresses  will    arise    if    thermal   expansion 
is    restrained.      The    stresses   in    a   laminate    are   the   result    cf 
thermal   and   load   stresses.      The    strength    d f    a    laminate    de- 
pends  upon    the   stress    state  and    the    failure    criterion    go- 
verning  failure. 

3 . 1  OBJECTIVES 

The  objectives  of  the  strength  investigation  are: 

1.  to  determine  the  stresses  in  a  laminate  due  to  any 
combination  cf  inpiane  loads,  bending  loads,  and 
"thermal  loads", 

2.  to  determine  the  strength  of  a  laminate  resulting 
from  these  stresses. 

3.2  DESCRIPTION    OF    THE    MODEL 

A    detailed   formulation    and    inscription    of    the    mathemati- 
cal   model    is    presented   in    a  technical    report     (8).      A    brief 
summary    follows. 

The   model    is    restricted    to    symmetric   laminates.      This   is 
hardly   a    limitation   on    the    model    because 

1.  nonsymmetric   laminates    are    extremely    uncommon,    and 

2.  symmetric   geometry    does   not    restrict    the    model    to 
symmetric   temperature   distributions. 

The   formulation    begins    with   the   formation    of    the    stiff- 
ness   matrix  of   a    lamina    in    the   natural   coordinates   of    the 
lamina.      The   stiffness    matrix    of    a   lamina-,?    is    obtained   from 
the    lamina    stiffness    matrices    by    coordinate    transformations 
from    the    natural   axes   of   a    lamina    to    the    laminate   coordinate 
system.      These   are   the    laminate    stress-strain    relations. 

The  strain-displacement    relations    (Kirohhof f-Love    theory) 
of    thin    plate   theory    are   combined    with   the    stress-strain    re- 
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laticns   to    give    the   laminate      force-displacement    relations. 
The    stresses    in    each    lamina  can    be   obtained    from    the    forces 
through   the   force-stress   relations.      Temperature    effects   en- 
ter   into    the   analysis   through   the    stress-strain    relations, 
which   in    turn   produced    "thermal    load"    terms    in    the   force- 
stress   relations. 

Once   the   lamina    stresses   are    known,    the    task    is   to    deter- 
mine  which,    if   any,    of   the    lamina    have    failed.      The   general- 
ized  Tsai-Wu   criterion    is    used   as   the   failure    criterion   for 
the    lamina.      The   failure   cf   one    or  more    lamina    ioes   not    im- 
ply   that    the    laminate    has    failed.      The   strength    analysis 
determines    the   ply    by   ply    sequence   of   failure.      After    a   ply 
(lamina)     has    failed,    the   ply    is    removed,    and   the    stress- 
strain  relations    for    the    modified   laminate    are    generated. 
This    process    is    repeated    until   all   lamina    have    failed.       If 
ail    cf   the    plies   have    failed,    the   laminate    has    failed.       If 
■*:he    loading   and   temperature   are    such    that    all    of    the    lamina 
have    not    failed,    then    the    laminate   has   not    failed.      In    this 
case    the    task   is   to   determine    "how   far"    the    laminate    is    from 
failure.       This   is    obtained   from    calculating   the    "loai    fac- 
tor". 

3.3  MODIFICATIONS    OF   THE    MODEL 

The   original    formulation   of   the   model    was    restricted   to 
symmetric    temperature    field  problems.       In    addition,    the 
analysis    was   restricted   to    problems   for    which    only   the    load 
factor  could    be    determined.      These   restrictions    have    been 
eliminated   in   the    present    model    cf   the    strength    analysis 
program    (STE5N)  . 

3.3.1        Npnsymmetric    Temperature    Fields 

The  strength    analysis   program     (STREN)     was    modified    to 
handle  nonsymmetric   temperature    distributions    which   result, 
for    example,    when   one   surface    cf    the    laminate    is    exposed    to 
a    fire.      Some   of   the    equations   presented    in    reference     (8) 
are    changed   as    fellows: 
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3.3.2        Ad^iiiSl^l  Failure   Factors 

In   it's    earlier   form,    STR2N      provided   the    load    factor    to 
failure    for   a    specified   temperature   distribution.      Tb.e    pro- 
gram   was    modified   to    handle   "temperature    factor",    and    "com- 
bined   factor"    problems    as    well.       each   of   these    problems    will 
be    briefly    described. 

In  the    load    factor    problem,    the    magnitude    of   temperature 
in    each    ply   is   specified   explicitly,    and   the    relative    magni- 
tudes  of    the    laminate   forces   and    bending    noments   are    given. 
If    the   laminate    load    vector    (Mx, Ny ,Nxy,dx, My, Mxy)    is    speci- 
fied  as    (1  ., .5, 0., 0. ,10.  ,0  .) ,    and  the   ioai    factor    is    deter- 
mined  to    be   R,    then    for    the  given   temperature    field,    +he    la- 
minate  will    fail    when    Nx=1.xR,    Ny=0.5xR,    Nxy=0.xR,    Mx=0.xR, 
My=10.xR,    and   Mxy=0.xR.      Reference    (8)    shows    that   the    load 
factcr   appears    as   the   scalar    multiplier   of    the    "lead    stress- 
es" . 

In   the    temperature    factor    problem,    the    laminate   load   mag- 
nitudes  are   specified   exactly,    and   only   the    relative    magni- 
tudes  of    the    ply   temperatures   are   known.       In    this   case,    the 
problem   is   to   determine   the  factor   which    when    multiplied    by 
the    relative    lamina   temperatures    gives   the    temperature    fieli 
at    laminate    failure. 

Finally,    we   consider   the  combined    factor    problem.       Here 
the    relative    maanitudes    cf    both    the    loads,    and   ply   tempera- 
tures  are    specified,    and   the   task   is    to    determine   the    "com- 
bined   factcr"    which,    upon    multiplication    with    each   of    the 
relative    loads   and    temperatures,    gives   the    loads   and   temper- 
atures  at    which    the    laminate    will    fail. 

In   each    case,    the    Tsai-Wu    failure    criterion    reduces    to 
the    guadratic   in   the    factor   R, 

aR2    +    3R   +    Y    =    0 

The    coefficients    a,   3 ,    and    Y  depend   upon    which    R    factor    is 
being    determined.      In    the    load    factor    case,  ex  ,3   ,    andY   ,    are 
the    a,    b,    and   c    of   Equations    (55)     through    (67)    of   reference 
(8) ,    respectively.      In    the    temperature    faotor   case,    a  is    b, 3 
is    a,    and    y  is   c.       In    the    combined    factcr    case, 
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a  =  a  +  b  -  Fx(aLx  +  aTx)  -  Fy(aLy  +  aTy) 
*  =  Fx(aLx  +  aTx}  "  Fy(aLy  +  ffTy> 


Y  =  -1. 


3.4  SAMPLE    M1MSES 

Reference    (8)    presents    some   results  of   the   strength    in- 
vestigation.     Those   results  showed  good   agreement    with    pre- 
vious  analyses   in   the    literature.      The   present    version    of 
STREN   was    used   to    generate    a    failure    locus    of    a    0/90    gra- 
phite-epoxy   laminate    at    a    temperature    of    200    F    below    the 
curing  temperature;    for    all   stress   states    of   the    form 
(Nx, Ny,Nxy,Mx, My, Mxy)  =  (Nx, Ny, 0,0 ,0,0,0,)  .       Since   there    are 
only    two    nor.    zero    stress    components    (Nx   and    Ny)  ,    the    locus 
of    failure    stress    states   can   be    shown   as    a    curve    in    tfx-Ny 
space.      This   curve   is   shown  in   Figure    15.       This    failure   lo- 
cus   was   generated    by    running   SIREN   for    37    combinations    of    Nx 
and    Ny    (    i.e.,     (Nx,Ny)     =    (1.,0.),     (1  .  ,  .5)  ,  .  .  .  ,  etc. )  . 
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Figure   27.    Failure   Locus    for    0/90    Graphite-Epoxy   Laminate 
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